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Abstract 

A magnetic field and the resulting Landau degeneracy enhance the infrared contributions to the 
quark mass gap. The gap does not grow arbitrarily, however, for models of asymptotic free 
interactions. For B — > oo, the magnetic field decouples from the dimensionally reduced self- 
consistent equations, so that the gap behaves as ~ Aqcd (or less), instead of ~ *J\eB\. On the 
other hand, the number of participants to the chiral condensate keeps increasing as ~ \eB\ so that 
1(^0)1 ~ |e-B|Aq CD . After the mass gap stops developing, nothing tempers the growth of screening 
effects as B — > oo. These features are utilized to interpret the reduction of critical temperatures 
for the chiral and deconfinement phase transitions at finite B, recently found on the lattice. The 
structures of mesons are analyzed and light mesons are identified. Applications for cold, dense 
quark matter are also briefly discussed. 



1. Introduction 

In past decades, systems in a magnetic field (B) have been useful laboratories to test theoretical 
ideas. A famous example is a system of cold atoms, in which a magnetic field controls the 
strength of the interactions. In QCD, similar utilities are also expected for the lattice Monte 
Carlo simulation at finite B [IH1]. In particular, we can study the nonperturbative gluon dynamics 
through polarization effects, controlling quark dynamics by a magnetic field. Such information 
may help the studies of cold, dense quark matter [5]. 

It seems that lattice studies already confirmed some of the theoretical ideas. At T = 0, a 
magnetic field enhances the size of the chiral condensate due to magnetic catalysis [HI [7] . A key 
feature of this phenomenon is the effective dimensional reduction. For B ^ 0, the phase space for 
the low energy particles and anti-particles is ~ |eS| J dpy, increasing the number of participants 
to the formation of the chiral condensate. This should be contrast to the B = case, where phase 
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space quickly decreases as ~ J \p\ 2 d\p\ in the infrared region. In this case, due to the small number 
of participants, the system needs sufficiently strong attractive forces to form chiral condensates. 

On the other hand, some surprises have been provided as well [5]. While B increases the 
chiral condensate below the (pseudo-) critical temperatures for the chiral restoration (T x ) and 
deconfinement (7b), those temperatures themselves decrease. This might contradict with our 
intuitions, if we think that a larger chiral condensate should generate a greater quark mass gap. 
Such thinking would suggest that (i) a larger quark mass gap should suppress thermal quark 
fluctuations, leading to increasing T x , and (ii) a larger quark mass gap suppresses quark loops, so 
that the results should approach to the pure gauge results, leading to increasing T D . 

To resolve this apparent contradiction, we shall argue that the quark mass gap at T = can 
stay around ~ Aqcd (or less) for large B. Then we can imagine that the gap at T < T X; d also 
stays around ~ Aq CD , without strongly suppressing thermal quark fluctuations and quark loops. 
If this is the case, the decreasing of critical temperatures would not be so unnatural. 

In addition, the aforementioned behavior of the quark mass gap does not contradict with 
the growing behavior of the chiral condensate, but instead naturally explains its -B-dependence 
at T = 0. In fact, the lattice results in |4j showed the behavior (^V')t=q ~ l e -B|A.QCD, for 
\eB\ > 0.3 GeV 2 ^> Aq CD (~ 0.04 GeV 2 ). Noting that the relation under the Landau quantization, 

(V>V)4D ~ \eB\ x (V^) 2 d , (1) 

we can see that {ipip)2D or the quark mass gap must be nearly ^-independent and O(Aqcd)- 

In this work we will carry out all the calculations in the large N c limitj^] The use of the 
large N c is motivated by at least three reasons: (i) At large N c , gluons are not screened, so the 
nonperturbative forces (i.e. the forces in the infrared) are stronger than the N c = 3 case. Such 
forces can be used to set the upper bound of the quark mass gap. (ii) The large N c limit has 
captured many qualitative aspects of the confined phase at B = 0. Therefore it is worth thinking 
and testing this approximation in the confined phase at finite B, since its validity and invalidity 
are not evident apriori. (iii) It is easy to imagine how the 1/N C corrections qualitatively modify the 
large N c results, and such corrections just provide welcomed effects for our scenario (see below). 

We will use the large N c limit to just claim that the quark mass gap does not grow much beyond 
Aqcd- To explain the reduction of the critical temperatures, in addition we have to argue the 
1/N C corrections. The quark loops as the 1/N C corrections screen the nonperturbative forces. As 
B increases, the screening effects become larger because more low energy particles can participate 
to the gluon polarization, due to the enhanced Landau degeneracy ~ \eB\ in the lowest Landau 



1 The large N c limit in a magnetic field was also studied in Ref. [9] from a different perspective from ours. 
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level (LLL) [8]. If the quark mass gap stops growing as suggested in our scenario, there is nothing 
to suppress the growth of the screening effects as B increases^} Therefore the nonperturbative 
forces are reduced at large B, and such reduction should lower the critical temperatures for given 
B. In addition, hadronic fluctuations as the 1/N C corrections also grow as B increases, helping 
the chiral symmetry to restore [10£J 

We will argue that the demanded (nearly) B-m dependent gap of O(Aqcd) can be derived, 
provided that it is dominantly created by the nonperturbative force mediated by the IR gluons. 
In particular, both the IR enhancement (that is more drastic than the perturbative 1/p 2 case) 
and the UV suppression of the gluon exchanges are crucial for our discussions. Since we will deal 
with the LLL which is essentially soft physics in the present paper, IR enhanced gluon is a key 
feature in this study (see [12] for a review). If we include only the perturbative 1/p 2 force, the gap 
is much smaller than Aqcd and depends on B at most logarithmically. Similar arguments have 
been used in studies of the quark mass function at finite quark density [131 [H] . 

To illustrate our points, we first consider the NJL model which does not have the abovemen- 
tioned properties. For \eB\ — > oo, the gap equation within the LLL approximation is 

M mL (B) = GtrS(x,x) — > G [ ^ - MnJl(jB) = A), (2) 

where f(q z ,B;A) is some UV regulator function. The contact interaction couples all states in the 
LLL so that the Landau degeneracy factor \eB\ for the LLL appears. The intrinsic property of 
the model is that the chiral condensate has the same 5-dependence as the mass gap: 

W>njl - -^M NJL (B). (3) 

Depending on the regularization schemes, Mnjl(-B) can be ~ |e_B| 1,/2 (proper time regularization 
[T5]). or ~ A (four momentum cutoff [ID]), or else. Each scheme has its own problems. For schemes 
predicting the growing behavior of the chiral condensate [16J, the quark mass gap also develops 
as B increases. Then at finite temperature, thermal quark contributions are largely reduced so 
that the increasing chiral restoration temperature is naturally expected. On the other hand, if the 
mass gap approaches to constant, the chiral condensate also does, contradicting with the lattice 
results. Therefore, as far as the relation like (|3| is retained, it seems that we have to abandon 
either the increasing chiral condensate or the reduction of critical temperatures. 

This dilemma can be bypassed if we use the gluon exchange type interactions with the IR 
enhancement and UV suppression. To emphasize the point, we consider a simple model for the 

2 This suggests that even at T = 0, the nonperturbative gluons will be screened out at some critical value of B 
such that the screening mass, mp ~ OsleBl 1 / 2 ~ N^ 1 ^ 2 \eB\ 1 / 2 , becomes comparable to Aqcd (see also SecjiJ. 
3 For hadronic fluctuations at small \eB\, see Ref. [TT], where chiral perturbation theory should be at work. 
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Figure 1: (a) The Schwinger-Dyson equation at large N c for the model in Eq.Q. (b) The distribution of states 
in the lowest Landau level for fixed p z . A state with momentum p± can strongly couple to states within a domain 
of |<7j_ -p±\ < A f ~ Aqcd- 

gluon exchange with these features (for the moment we ignore spinor structures), 

D(q)=G6(A 2 f -q 2 ), (A f ~ A QCD ) (4) 

which was proposed in [14J. In this model, the quark mass function appears to be momentum 
dependent. For B — > oo, the Schwinger-Dyson equation at large N c (Fig Ilk) is 



M(p; B) ~ G J 0^ tr S% L (q z ) 6{A) - \p- q\ 2 ) , (5) 

where B is used to reduce the quark propagator to the (l+l)-dimensional one and to separate 
higher Landau levels from the LLL. In contrast to the previous case, the factor \eB\ does not 
appear in front of the integral. This is because the interaction (|4]) does not couple all the states 
in the LLL, but couples the states having similar momenta (Figjlb). This feature makes the gap 
B- independent. In fact, carrying out the integral over the transverse momenta, we get 



M(p; B)*£[ ^ 1 m m ~ \P* ~ 1*?)F{p z - g.) , (0) 

2tt J In y/qi + M 2 (g; B) 



where F(k z ) = Aj — \k z \ 2 . Note that the equation does not have any explicit B dependence, so 
the mass gap is solely determined by the scale A^, i.e. M(p; B) = M\ f (p). 

Another important feature is the damping behavior of the mass function at large p. To see 
this, in ^ we take p z — > oo, for which q z must go to oo. Then the integrand in the RHS goes to 
zero. The phase space integral cannot compensate the damping behavior of the integrand, because 
the integral is limited within the finite domain A^ around p z . So M(p z ) — > as p z — > oo. Thus 
there must be some damping region of M(p z ) which we will denote A x as a function of Af. The 
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emergence of the damping scale A x makes the chiral condensate UV finite^ 

<^> ~ - \eB\ ^ - M(9z) ~ - \ e B\ F(A f ) , (7) 
J-a x v^F + ^fe) 

where F(Af) is a function of Ay. The chiral condensate depends on B linearly, as advertised. 

Note that all of the discussions presented so far did not ask whether the interaction is confining 
or not. All we needed was the IR enhancement. Nevertheless, it is interesting to further investigate 
a model which not only contains properties of the IR enhancement, but also captures certain 
aspects of the confinement. For this purpose, we will use the model of Gribov-Zwanziger [17] 
which contains the linear rising potential between colored charges. Confinement is expressed as 
the absence of the quark continuum in the meson spectra. We will analyze the self-consistent 
equations for the quark self-energy and meson states at large N c and large B, by reducing them 
into those of the 't Hooft model. Then we will rederive the aforementioned conclusions. 

2. Dimensional reduction 

We consider the Euclidean action (convention: g^ v = 5^, 7^ = 7^): 

S E = J d 4 x $ [0 + iQA + m\ ij + S int , (8) 

where A is a U(l) em gauge field, Q is a flavor matrix for electric charges, and m is the current 
quark mass matrix. The color gauge interaction is treated as S int . 

We apply an external, uniform magnetic field in spatial 3-direction^J which can be given by a 
vector potential (Ai, • • • , A4) = (0, Bxi, 0, 0). For later convenience, it is useful to introduce the 
projection matrices. We decompose fermion fields, 

PS = 1 ± Sgn !f B)(T3 , ^ = 17172, V± = ^, a^± = ±sgn(QB)^ ± . (9) 
We expand the fermion fields in bases which diagonalize the unperturbed Lagrangian, 

Mx) = J2 I ^Ml,P)e~ m H l (xi,P2), (10) 



where p = (0,p 2 ,p 3 ,p 4 ). Hi(xi;p 2 ) is the harmonic oscillator base with mui = \QB\, whose center 
is located at X\ = pijQB. If we want to write the Lagrangian as a sum of the Landau levels, we 



4 In vacuum, A x effectively plays the role of the UV cutoff in the NJL model. 



5 We try to explain only the lattice results for uniform B in the quenched QED. If we had some charged 
condensates (see recent discussions, [18— 20J), the dynamical QED would generate vortices in electric superconductor. 



5 



should relabel tp+(n, x) — > xp n+ (x) , ip_{n — 1, x) — )■ ip n -(x) , t/>+(0, x) — > . Now the index 

n characterizes the Landau level. The unperturbed action becomes 



/ d 3 X lc (x) + Y,£n(l 
J \ n=l 



d 4 X Ainpert. (x) = I d*X \ Cq(x) + ^> " CJx) \ , (11) 



after carrying out the integral over x\. The Lagrangian density is (7^ = (73, 74) and 7j_ = (71, 72)) 
£o = x(x)(h + ™)P+x(x)i £n = Mx) ($ L + isgn(QB)y/2n\QB\ l2 + m) if> n (x) . (12) 
The propagators for Landau levels are 



x ' J (2it) A p L + 1m 

We also expand the color gauge interactions in the Ritus bases instead of usual Fourier bases [2"T| . 

d 4 x ihptatiAfnix) = 22 / (2 , 6 i>(l,P)lnt a i/;(l',q)A l l ia (p-q;p2,q2), (14) 

where the gauge field is convoluted with the harmonic oscillator bases, 

4jl(p-g;P2,g 2 ) = J ^ A^p - q,h) J dx x H l {x x -p 2 )H v {x x -q 2 )e- ik ^. (15) 

For later convenience, it is useful to prepare the propagator for this convoluted form, 

A l laiP -T,P2,q2) A 3 2 (q - p; q 2 ,p 2 
= / j dx x j dy 1 H l {x 1 -p 2 )H l ix 1 -q 2 )Df v {p-qM)H (16) 

where we used (A^ a (p)A ub (q)) = (2n) 4 5 4 (p + q)D a ^ y (p). Let us note that the Hi(xi,p 2 ) depends 
on coordinates in the combination, x\ — p 2 /QB. We can simultaneously shift x x — p 2 /QB — > xi, 
and yi — p 2 /QB — > y x , and get 

Hi(x 1 ;p 2 )Hi>(x 1 ;q2)Hj / (y 1 ;q 2 )H j (y 1 ;p 2 ) ->■ H^x^Hi^Xf, q 2 -p 2 )H f (yr, qi -p 2 )Hj{yx) , (17) 

6 It is convenient to attach the Ritus wavefunctions to the vertices instead of propagators, because the different 
Landau orbits couple only after interactions are turned on. 



keeping the other part in Eq.(16) invariant. Thus the integral depends on p 2 and q 2 only through 
Vi ~ <?2- This is natural consequence of the gauge invariance, A 2 = Bx\ — > B(x\ — c), that affects 
the origin of p 2 and q 2 but should not affect the final result. 

This complicated expression can be drastically simplified if (i) the interaction shows an IR 
enhancement around ~ Aqcd, while it damps quickly in the UV, and (ii) the magnetic field is 
sufficiently strong, \QB\ 3> Aqcd- The gluon propagator rapidly damps for p — q and k\ much 
larger than Aqcd, cutting off the domain of integral. Note also that the function H(xi) contains 
the Gaussian function, exp[— \QB\x\/2], cutting off the domain of X\ by ~ 1/ y/\QB\, as far as we 
consider not very high Landau orbits. 

Assembling all these properties, we can make replacements at very large B, 

e - ikl ( Xl - yi) ^ 1; (|A:a(ar 1 -y 1 )|~AQCD/v1W<l) 
Hi(x 1 ;p 2 - q 2 )H v (xi) -> Hi(xi)H v (xi) , [\p 2 - q 2 \ ~ A QC d < V\Q B \ 



which simplify the Eq.(16) as 



A^ a (p-q;p 2 ,q 2 )A^(q-p;q 2 ,p 2 )) — ► J ^ S lv S fj D%(p - q, h) . (19) 

(Note that indices are for the orbital quanta. Whether they coincide with the Landau level indices 
depends on 7-matrices to which gluons couple, because 7^'s flip spins.) 

The final expression itself coincides with a naive expectation and is not surprising. Ultimately, 
as B — > 00 a hopping from one Landau orbit to others (transverse dynamics) is completely sup- 
pressed. A nontrivial consequence of the asymptotic free theories, however, is that the separation 
of the Landau levels is achieved at a relatively small magnetic field, compared to forces like 1/p 2 . 
For the 1/p 2 force, the damping of the UV force is much milder so that the LLL and higher 
levels do not decouple quickly. The validity of separation is quantified by examining how the 



replacements in Eq.(18) work. 



3. The self-consistent equations at large iV c within a confining model 

Now we consider the self-consistent equations, the Schwinger-Dyson and Bethe-Salpeter equa- 
tions. Below we set m = 0. We use the gluon propagator of the Gribov-Zwanziger type |17j . 



D; b u (k) = -^D, u (k), J D 44 (fc) = -5 44 Ay(*)=iVAO = 0. (20) 

which is motivated by Coloumb gauge studies. D M gives a linear rising potential for the color 
charges, and we dropped off terms without the strong IR enhancement, a is a string tension of 
0(Aq CD ), Cf = (N 2 — 1)/2N C is Casimir for the adjoint representation. 
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We will assume that the self-energies are diagonal for each Landau level, and are translational 



invariant. Then the Schwinger-Dyson equation, after applying our approximations (19), becomes 

d 3 g 



$l + £_l + T, m )(n,p) = J -^4^ i7 4 5'(ra,g;S)i74 J^D 44 (p-q,k 



(21) 



that is diagonal for the Landau orbits. 

The computations are particularly simple for the LLL. Below we drop the subscript for the 
LLL. The matrix 7_|_ drops off because of the projection operator P®, so that RHS of the equation 
does not have 7j_, meaning that = in LHS. Then the equation looks like 

d 3 g i[4 L -$ L (g)] +S m (g) ^ fdh 
[ gi -S L (g)] 2 + S^(g) 



+ ^ ® = / rda ^ r:;^^z 174 *r D ^ p - *• h) • (22) 



Recall that the quark transverse momentum is gauge dependent in the sense that A 2 = Bx\ — > 
B(xi — c) affects the origin of the transverse momenta as p 2 —> p 2 + QBc and q 2 —> q 2 + QBc. Doing 
this shift, the gluon propagator is unaffected while the modification appears only through the self- 
energy as £(p£,p 2 ) — > S(pi,p 2 + QBc), and £(pz,, q 2 ) — > £(pl, q 2 + QBc). This replacement does 
not change the structure of the self-consistent equation at all, so that we obtain the same solutions 
for the two cases, £(??l> P2) = ^(pl,P2 + QBc), which means that E(p) is p 2 independent. 
As a consequence, we can factorize the equation: 

/ v ±v v x f d 2 q L . i[<j L -$ L (q L )] +S m (g L ) , f dhdq 2 

(tL + Em)M = J W? 174 171 J W " ? ' fe) ■ (23) 



The gluon propagator is smeared by integrating out the transverse momentum, 

dkidq 2 .„ _ f d 2 k± —871a 2a 



which is a confining propagator in (1+1) dimensions. The form of the Schwinger-Dyson equation 
is exactly same as that of 't Hooft model |22J in axial gauge, A z = 0, whose solution is known. 
The string tension is related to the two dimensional gauge coupling as N c g 2D = 4a. 

The magnetic field disappears from the equation. The only role of the magnetic field is to 
make the unperturbed quark propagator (l+l)-dimensional and to separate the LLL from the 
other Landau orbits. The only scale in the equation is a ~ Aq CD . Using the dressed (1+1)- 
dimensional quark propagator for Eq.([7]), we can get the chiral condensate at finite B. 

Since the quark self-energies are not always well-defined, we will also consider the meson states 
which are free from any ambiguities. One can estimate an effective or constituent quark mass by 
examining the meson spectra. To do this, we study the Bethe-Salpeter equation. 
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We shall move to Minkowski space, and treat the homogenous Bethe-Salpeter equation assum- 
ing that the total momentum of a quark and an anti-quark is sufficiently close to the pole. The 
equation for a color singlet channel is 

*{P\ *)% = - J [ S(k - P) 70 A;) 70 S(k + P) } f ^D m {k - q) , (25) 

where \P(P; q) is a meson wavefunction with total momentum P and relative momentum q for a 
quark and an anti-quark. The indices a, (3 (/, /') are for spinor (flavor) indices. We can decompose 
the wavefunction into different spinor combinations and flavor structures, 

= (*s + *57s + ^l7l + • • • ) ap Iff + ■■■ ■ (26) 

We must apply further decompositions or take appropriate linear combinations, because a mag- 
netic field breaks isospin and rotational symmetries explicitly. Below we restrict ourselves to the 
meson states with electric charges and spinor combinations for which we can close the equa- 
tion only by the LLL. One can select out such mesons using the projection operator P+ = 
(1 + isgn(Q_B)7i7 2 )/2. For the two flavor case, the following structures satisfy the condition, 



(uu,dd)®(l,j 5 ,j L ,j L j 5 ,a L L',cr±±>), (ud,du)®(j±,j±j 5 ,a L ±). (27) 

For instance, neutral pions remain light, while charged ones acquire the masses ~ y/\eB\. Other 
examples are vector mesons. The longitudinal (transverse) component of neutral (charged) vector 
mesons can remain light, while others nolQ This observation seems to be consistent with recent 
lattice results [19] and model calculations [23] . 

Below we consider the neutral scalar component, as an illustration, taking only the LLL 
into account. As we saw, the dressed quark propagator is independent of the transverse momen- 
tum. Then we can conclude that *&s(P] o) is independent of P±, because Pj_-dependence appears 
only through so that different P± gives the same equation and thereby the same solution. We 
can also conclude that ^s(P',q) is independent of q±, because simultaneous shift of momenta, 
q± — > q± + c and k± k± + c, does not affect the equation. 

Therefore we can again factorize the equation, 

*5(^x;?l) = - J ^y 2 [ S ( k L- p L)lo^s(PL;k L ) lo S(k L + P L )} J ^D 00 (k-q), (28) 



7 More careful considerations are necessary beyond the large N c limit where the annihilation diagrams may 
contribute. In this respect studies of instantons in a magnetic field are important [23] , 
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QCD 




Figure 2: A structure of mesons listed in Eq.(27). With a magnetic field, a quark and an anti-quark aligns along 



direction of the magnetic field. The meson can move in the transverse direction without costing much energy. 



and obtain the Bethe-Salpeter equation in the 't Hooft model. There is no B-dependence in the 
equation, so masses of mesons are characterized by current quark masses and Aqcd- There are an 
infinite tower of meson spectra, due to the string oscillations along ^-direction. 

The total momentum of a quark and an anti-quark is independent of transverse momentum, 
meaning that the meson can move freely in the transverse direction without costing energ}]^} The 
inclusion of the higher Landau level can generate the dependence on the transverse momentum, 
which is suppressed by a factor of Aqcd /lei?! 1 / 2 . This leads to a spatially anisotropic meson 
Lagrangian [8] for which mesonic fluctuations are strong due to enhanced phase space in the 
infrared. It may significantly affect dynamics at finite temperature |25| . like a system close to the 
quasi-long range ordeiQ [27] . 

On the other hand, the independence of the relative momentum means that the internal meson 
wavefunction behaves as 



(29) 



that is, a quark and an anti-quark align along the ^-direction, and the inter-particle distance 
is ~ 1/Aqcd (Fig|3]) . Couplings with the higher Landau levels will introduce the width for the 
transverse wavefunction, and the mean square radius should be a/ (fj_) ~ Aq CD /|Q-B|. 



4. Conclusion 

We have discussed that the quark mass gap in the presence of a magnetic field stays around 
~ Aqcd or less, provided that it is mainly generated by the nonperturbative part of the gluon 
exchange. Accordingly the chiral condensate grows at most linearly as a function of B. This 



8 The meson dynamics is effectively (l+l)-dimensional, in the sense that the transverse momenta do not affect 
much the energy dispersion. This means that mesons can easily propagate in the transverse directions. 

9 In the strict quasi-long range order, ("ipip) ~ p(e 10 ) — > due to the IR divergence of the phase fluctuations. 
But the quark mass gap can remain finite as far as p ^ [26| . 
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Figure 4: The momentum space structure near the Fermi surface. If the curvature is small, the structure of the 
Schwinger-Dyson equation for a non-uniform chiral condensate is very similar to that for the LLL in a magnetized 
system (See also Fig[lb). 

tendency seems to match with the current lattice data. But to complete our arguments, more 
detailed discussions about the UV tail of the gluon propagator are necessary. We leave it for 
future studies. 

We have not explicitly taken into account the screening effects for gluons by taking large 
N c . Including 1/N C corrections, the screening effects should grow as B increases, reducing the 
nonperturbative forces. Accordingly the growth rate of the chiral condensate should become 
smaller than the large N c estimate. We expect that for an extremely strong magnetic field, the 
screening effects will exceed some critical strength to make the system deconfined (Fig|3]). For 
T = 0, the critical B is roughly estimated as |ei?| 1//2 ~ A^" c 1//2 Aqcd, by equating the numbers of 
virtual quark and gluon excitations at low energies, iV c |e-B|AQCD ~ -^c-^-qcd- 

This situation seems to be quite analogous to what would happen in cold, dense quark matter. 
At large quark chemical potential fi q , the area of the quark Fermi surface is enhanced. Accordingly 
the phase space for quark excitations at low energy increases as ~ Nc/i^Aqcd, so that the screening 
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effects for gluons become strongei^ The critical chemical potential for deconfinement has been 

1 /2 

estimated to be /i g ~ Av Aq CD [2H]- Detailed quantitative estimates crucially depend upon the 
existence of the quark mass gap near the Fermi surface. The gap of O(Aqcd) can emerge if we take 
into account the non-uniform chiral condensates. The corresponding Schwinger-Dyson equation 
[H] takes very similar form as that in a magnetized system, Eq. ^ . This correspondence becomes 
better at larger fj, q , for which the curvature of the Fermi surface is small enough (Figji]). 

In summary, the lattice Monte Carlo simulation at finite B is a promising tool to extract out 
quantitative information about the gluon polarizations. Such information has great relevance to 
understanding dense quark matter, or more concretely, physics near the Fermi surface. On the 
other hand, information about the bulk Fermi sea should be supplemented by the lattice studies 
for the two-color QCD (or isospin QCD) at large quark (isospin) density without a magnetic field 
j. These discussions will be expanded elsewhere. 
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